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Abstract. Rate coefficients for inelastic collisions between Li and H atoms covering all transitions between the asymptotic 
states Li(2s,2p,3s,3p,3d,4s,4p,4d,4/)+H(ls) and Li~''+H' are presented for the temperature range 2000-8000 K based on 
recent cross-section calculations. The data are of sufficient completeness for non-LTE modelling of the Li I 670.8 nm and 
610.4 nm features in late-type stellar atmospheres. Non-LTE radiative transfer calculations in both ID and 3D model atmo- 
spheres have been carried out for test cases of particular interest. Our detailed calculations show that the classical modified 
Drawin-formula for coUisional excitation and de-excitation (Li*+H^Li*'+H) over-estimates the cross-sections by typically 
several orders of magnitude and consequently that these reactions are negligible for the line formation process. However, the 
charge transfer reactions coUisional ion-pair production and mutual neutralization (Li*+H^Li"'"+H^) are of importance in 
thermalizing Li. In particular, 3D non-LTE calculations of the Li I 670.8 nm line in metal-poor halo stars suggest that ID 
non-LTE results over-estimate the Li abundance by up to about 0. 1 dex, aggrevating the discrepancy between the observed Li 
abundances and the primordial Li abundance as inferred by the WMAP analysis of the cosmic microwave background. 
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1. Introduction 

Lithium abundances in stellar atmospheres are key observa- 
tional parameters in astrophysics. They tell us much about stel- 
lar evolution and nucleosynthesis, Big Bang nucleosynthesis, 
and the role of cosmic rays in Galactic chemical evolution, to 
name a few (see Lambert ll9931 Carlsson et al. 119941 and refer- 
ences therein). Stellar Li abundances are in practice determined 
from the 670.8 nm resonance line of Li I and in a few cases the 
610.4nm subordinate line. Thus, the ability to accurately inter- 
pret these lines is important. As Li I has a low ionization poten- 
tial, non-local thermodynamic equilibrium (non-LTE) effects 
can be expected to be important (Steenbock & Holweger[1984| 
Carlsson et al. .l994. Asplund et al. 2003 1. 

Steenbock and Holweger ([p84i were the first to suggest 
that inelastic collisions with hydrogen might be important in 
thermalizing Li L In photospheres of late-type stars neutral hy- 
drogen atoms outnumber electrons by typically TVjy /A^e ~ 10*, 
even more in metal poor stars. Thus, weight of numbers may 
overcome the expectation that electron collisions are a more 
efficient thermalizing mechanism firstly due to their higher 
thermal velocity and therefore higher collision rate, and sec- 
ondly since for speeds of interest electron cross sections are 
generally larger than those for neutrals since the neutral col- 
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lisions are adiabatic while the electron collisions are not (e.g. 
Lambert 1993'). In such photospheres kT ^ 0.2-0.6 eV and 
thus for typical optical transitions the low energy collisions 
just above the threshold are most important in determining 
the collision rates. Steenbock & Holweger, in the absence of 
any better experimental or theoretical guidance, calculated the 
coUisional excitation and ionization rates by modifying a for- 
mula from Drawin f 1968i for H+H which is itself a modi- 
fication of the classical Thomson formula for excitation by 
electrons (cf. Lambert [1993}. In recent years some work has 
been devoted to improve these highly uncertain H collision 
cross-sections. Low-energy experimental data have been ob- 
tained (Fleck et al. ll991> and ab initio calculations performed 
down to the threshold (Belyaev et al. "1999^ for the Na-nH 
system with good agreement between theory and experiment. 
However, these results are for only the very lowest states of this 
one system. For non-LTE calculations, estimates are required 
for transitions between all states which might affect the popu- 
lation of the states of interest. It is for this reason that the mod- 
ified Drawin formula (often with a scaling factor determined 
astrophysically) is still often in use despite that it is known to 
overestimate the Na(3s)H-H— >Na(3p)H-H collision rate by sev- 
eral orders of magnitude. 

Recently, calculations for Li-nH collisions for all transi- 
tions between the states Li(2s,2p,3s,3p,3d,4s,4p,4(i,4/)H-H(ls) 
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and Li++H^ have been performed (Belyaev & Barklem!2003^. 
This paper should be consulted for details and justification of 
the approaches used. In brief summary, on the basis of the inter- 
pretation of the experimental and theoretical results for Na+H 
(Fleck et al. ll991l Belyaev et al. ll999> . non-adiabatic regions 
associated with avoided ionic crossings in the adiabatic poten- 
tial curves are expected to provide the basic mechanism for the 
transitions at low energy. A simple physical interpretation of 
this process is that at certain separations during the collision the 
optical electron associated with the Li atom may tunnel to the 
H atom, leading to a predominantly ionic charge distribution 
(Li+H^) of the LiH quasi-molecule. Later, at another avoided 
crossing, the electron may tunnel back towards the Li atom 
into a different covalent charge distribution leading asymptot- 
ically to a different final Li state. The electron may also re- 
main with the H atom leading to an ion-pair production reac- 
tion (Li*H-H— >Li+H-H~). For the transitions between the lowest 
four states a fully quantum mechanical approach has been used, 
while for the transitions involving higher states, a multichannel 
Landau-Zener model has been used. The model approach is 
justified for the higher crossings where the requirements of the 
Landau-Zener model are better fulfilled. 

For the charge transfer reactions Li*H-H^Li+H-H^, re- 
ferred to as ion-pair production and mutual neutralization, our 
calculations use the model approach. More elaborate calcula- 
tions, using a diabatic molecular expansion method, have been 
made by Croft et al. ( 1999b i. This process has a remarkably 
large cross section near the threshold in aUcali-hydrogen sys- 
tems, particularly to and from the first excited s state. This 
has been predicted theoretically (Bates & Boyd 1956 Janev & 
Radulovic .l978i Croft et al. 1999a. Dickinson et al..l99 9 1 and 
is consistent with experiment (Peart & Havton [l994> . Despite 
that this process has been known to be remarkably efficient in 
various fields (e.g. plasma physics, early universe chemistry) 
for some time, to our knowledge this is the first time it has 
been considered in a stellar atmosphere application of Li. 

2. Collision rate coefficients 

The rate coefficients (av) for excitation and de-excitation 
H(ls)H-Li(nO— *H(ls)H-Li(n7'), and the charge transfer reac- 
tions Li(n/)H-H^Li+H-H~ are presented in Tabled The data 
have been obtained by integrating a Maxwellian velocity dis- 
tribution with cross-sections from Belyaev & Barklem (2003) 
in the case of excitation and de-excitation, and from Croft 
et al. ( fT999b.) in the case of mutual neutralization, with ion-pair 
production cross sections obtained via the principle of detailed 
balance. Thus the data represent, in our opinion, the best avail- 
able estimates. Note that all forward and backward rates were 
computed separately, but should be connected by the detailed 
balance relation. Inspection of the data reveals the remarkably 
large rate coefficients for Li(n/)H-H^Li+H-H^, particularly for 
the n = 3 states. 

For the temperatures here the uncertainty in the rate coef- 
ficients is determined by the uncertainty in the cross sections 
near the threshold. The uncertainty in the input cross sections 
is discussed in detail in the relevant papers. In summary, the ac- 
curacy for transitions between the lower states (2s, 2p and 3s) is 



estimated to be better than a factor of 2. For the excitations in- 
volving higher states the data are expected to be accurate within 
an order of magnitude. Comparison with experimental results 
for the total neutralization cross-section with D^ (which uses 
the same interaction potentials as for H^) finds the theoretical 
result about 20% greater at the lowest measured energy 0.68 
eV (see Croft et al. I1999al l. with a general trend of the dis- 
agreement becoming larger at lower energies, recalling that the 
collision energies of importance here are E w 0.2 — 0.6 eV. 

Lambert J1993t presented unpublished estimated rate co- 
efficients for the 2s^2p transition computed by Allen & 
Dickinson. Our data are roughly an order of magnitude smaller. 

Comparisons were made with the modified Drawin formula 
for all optically allowed transitions at 5000 K. Except for an ex- 
treme outlier (3s^4p), the Drawin formula gives a rate coeffi- 
cient between roughly one and six orders of magnitude greater 
than our result. This reinforces the view that the Drawin for- 
mula typically greatly over-estimates the collisional efficiency, 
as judged by the available (admittedly rather scant) evidence. 

3. Application to non-LTE formation of Li lines 

To assess the impact of the data on the interpretation of 
stellar spectra, we have performed non-LTE calculations for 
Li in cases of particular interest using the MULTI-code for 
ID statistical equilibrium problems (Carlsson 1986) and the 
MULTl3D-code (Botnen 1997; Botnen & Carlsson 119991 
Asplund et al..2003) for the 3D case. The adopted model atom 
is the 21-level atom compiled by Carlsson et al. J1994> . To 
test the impact of the new data three model atoms were used, 
namely the original model atom with no H collisions (noH), 
the model with the new H collision data added (wH), and a 
model with the H collision rates halved (wO . 5H); the LTE re- 
sults have been obtained with the same codes using an atom 
with extremely large collisional cross-sections to ensure con- 
sistency with the non-LTE calculations. The ID model atmo- 
spheres are from the MARCS code (Gustafsson et al. 1975 and 
subsequent updates). The 3D models and calculations are as 
detailed in Asplund et al. (1999' '2003). Calculations were per- 
formed for the Sun, the metal-poor subgiant HD 140283 and 
the metal-poor turnoff star HD84937. Equivalent widths for the 
670.8 nm line are presented in Tabled Departure coefficients 
for the lower states in the ID non-LTE models for the Sun and 
HD140283 are presented in Fig.[Tl 

Although quantitatively different, the overall effects of the 
H collisions are similar in the ID and 3D cases. Almost all of 
the difference resulting from the inclusion of the H collision 
data is caused by the reaction Li(3s)-)-H^Li"*"-)-H~. Removing 
all H collisions but this reaction results in equivalent widths 
within 1% of the results for the wH model. This reaction 
causes a stronger collisional coupling between the ionization 
stages, nearly thermalizing the 3s state. Photon suction (Bruls 
et al.[l992j results in the other levels being affected, in partic- 
ular the ground state is pushed towards LTE, or even into over- 
population as seen in the case of HD140283 in ID, at around 
log T500 = —2. As clear from Table|2] the results are not partic- 
ularly sensitive to the details of the H collisional cross-sections 
with only small differences resulting when decreasing all H 



p. S. Barklem et al.: Inelastic H+Li and H +Li collisions in stellar atmospheres 



3 



Table 1. Rate coefficients {av), in units of cm'^ s ^, for the processes H(ls)+Li(7iO— >H(ls)+Li(n'Z'), or where indicated 
H(ls)+Li(nO^H-+Li+(ls2) and H-+Li+(ls2)->H(ls)+Li(n'/'), for selected temperatures in the range T = 2000-8000 K. 
Asterisks indicate transitions where the calculated cross sections are considered unreliable. 



Initial final state n l 



state nl 2s 2p 3s 3p 3d 4s 4p 4d 4f H~+Li 

2000K 





2s 








9 


18E — 


20 


1 


20E — 


25 


1 


67E— 


27 


4 


16E — 28 


2 


18E — 


29 


2 


19E— 


30 


1 


85E— 


30 


1 


84E— 


30 


2 


HE — 


25 




2p 


1 


39E— 


15 








8 


31E — 


18 


2 


76E— 


19 


7 


91E — 20 


6 


23E — 


21 


7 


04E— 


22 


4 


66E— 


22 


3 


26E— 


22 


2 


06E— 


19 














74E^ 


13 


















































3p 


2 


56E- 


18 


2 


80E- 


14 


9 


78E- 


11 








2 


23E-10 


3 


55E- 


11 


3 


30E- 


12 


2 


21E- 


12 


2 


27E- 


13 


4 


25E- 


11 




3d 


4 


91E— 


19 


g 


22E— 


15 


\ 


42E— 




2 


72E— 


10 






2 


92E— 


12 


2 


74E— 


13 


2 


83E— 


13 


2 


63E— 


14 


3 


27E— 


12 




4s 


\ 


TOE— 


18 


3 


53E— 


14 


5 


15E— 


10 


2 


97E— 


09 


2 


lOE— 10 


























3 


42E— 


13 






\ 


7SE— 


19 


2 


97E— 


15 


4 


51E— 




2 


58E— 


10 


2 


72E— 11 
















* 






J, 




3 


95E— 


14 




4d 


9 


07E— 


20 


\ 


30E— 


15 


\ 


89E— 




g 


67E— 




7 


31E— 12 


































4f 


2 


ggj; 


20 




gyg 


16 


3 


^QJ^ 


12 




49£ 


12 




2Qjr 23 
































H 


+L,i ' 


4 


12E — 


13 


2 


65E— 


11 


1 


07E— 


07 


5 


39E— 


08 


5 


34E— 09 

4t)UUK 


7 


63E— 


12 


9 


30E— 


13 






















2s 








3 


62E — 


17 


2 


70E — 


21 


9 


94E— 


23 


1 


58E — 23 


2 


30E— 


24 


3 


88E— 


25 


3 


46E — 


25 


3 


45E— 


25 


4 


62E — 


19 




2p 


2 


57E— 


15 








1 


23E — 


15 


7 


06E— 


17 


1 


34E— 17 


3 


lOE — 


18 


5 


86E— 


19 


4 


09E— 


19 


2 


87E— 


19 


7 


92E — 


16 










17 


3 


09E — 


13 








5 


95E— 




9 












79E— 












66E— 


13 


5 


53E — 


10 




3p 


2 


25E- 


18 


2 


25E- 


14 


7 


56E- 


11 








1 


26E-10 


5 


96E- 


11 


8 


73E- 


12 


6 


14E- 


12 


6 


36E- 


13 


3 


42E- 


10 




3d 


2 


42E— 


19 


2 


91E— 


15 


7 


93E— 


12 


g 


56E— 


22 






4 


35E— 


12 


g 


38E— 


13 


4 


49E— 


13 


4 


Q2E— 


14 


2 


39E— 


22 




4s 


g 


07E— 


19 


]^ 


27E— 


14 


]^ 


92E— 


10 


7 


62E— 


10 




16E— 11 


























2 


47E— 


12 




tp 


g 


33E— 


20 


\ 


06E— 


15 


\ 


61E— 


22 


5 


69E— 


22 


g 


19E— 12 






















J, 




2 


54E— 


13 




4d 


3 


23E— 


20 


4 


63E— 


16 


g 


75E— 


12 


2 


39E— 


22 


2 


62E— 12 


































4f 


9 


60E— 


21 


2 


09E— 


16 


]^ 


1 IE— 


12 


2 


98E— 


12 


2 


84E— 13 
































H 


; + 
+Ll 


1 


29E — 


12 


3 


HE — 


11 


8 


67E— 


08 


4 


22E — 


08 


4 


oob— uy 

6000K 


1 


39E— 


11 


2 


60E— 


12 






















2s 








4 


18E — 


16 


1 


17E— 


19 


5 


33E— 


21 


5 


9Uil — 22 


8 


37E— 


23 


1 


68E— 


23 


1 


52E — 


23 


1 


52E— 


23 


8 


OSE— 


17 




2p 


4 


97E — 


15 








9 


51E — 


15 


5 


53E— 


16 


7 


32E — 17 


1 


89E— 


17 


4 


24E— 


18 


3 


02E — 


18 


2 


12E— 


18 


1 


55E— 


14 














45E^ 


13 








* 


30E- 








































3p 


2 


96E- 


18 


2 


58E- 


14 


6 


75E- 


11 








8 


37E-11 


5 


46E- 


11 


9 


29E- 


12 


6 


64E- 


12 


6 


91E- 


13 


6 


59E- 


10 




3d 


2 


13E- 


19 


2 


23E- 


15 


5 


74E- 


12 


5 


44E- 


11 






3 


82E- 


12 


6 


51E- 


13 


4 


65E- 


13 


4 


18E- 


14 


4 


67E- 


11 




4s 


3 


32E- 


19 


6 


94E- 


15 


1 


06E- 


10 


4 


27E- 


10 


4 


59E-11 


























3 


22E- 


12 




4p 


3 


46E- 


20 


5 


76E- 


16 


8 


81E- 


12 


3 


12E- 


11 


3 


39E-12 


























3 


36E- 


13 




4d 


1 


76E- 


20 


2 


53E- 


16 


3 


68E- 


12 


1 


31E- 


11 


1 


44E-12 


































4f 


5 


24E- 


21 


1 


14E- 


16 


6 


04E- 


13 


1 


08E- 


12 


1 


OlE-13 
































H 


-+Li+ 


2 


54E- 


12 


4 


llE- 


11 


7 


95E- 


08 


3 


74E- 


08 


4 


05E-09 
8000K 


2 


29E- 


11 


5 


05E- 


12 






















2s 








1 


90E- 


15 


1 


llE- 


18 


5 


49E- 


20 


4 


86E-21 


4 


44E- 


22 


9 


73E- 


23 


8 


91E- 


23 


8 


89E- 


23 


1 


16E- 


15 




2p 


9 


27E- 


15 








3 


24E- 


14 


1 


84E- 


15 


1 


92E-16 


4 


12E- 


17 


1 


OOE- 


17 


7 


22E- 


18 


5 


06E- 


18 


7 


81E- 


14 




3s 


1 


48E- 


16 


8 


89E- 


13 








9 


90E- 


11 


1 


13E-11 


1 


71E- 


11 


3 


36E- 


12 


2 


34E- 


12 


5 


14E- 


13 


3 


04E- 


09 




3p 


4 


77E- 


18 


3 


28E- 


14 


6 


44E- 


11 








6 


14E-11 


4 


60E- 


11 


8 


44E- 


12 


6 


08E- 


12 


6 


34E- 


13 


9 


02E- 


10 




3d 


2 


69E- 


19 


2 


19E- 


15 


4 


73E- 


12 


3 


91E- 


11 






3 


16E- 


12 


5 


79E- 


13 


4 


17E- 


13 


3 


76E- 


14 


6 


62E- 


11 




4s 


2 


16E- 


19 


4 


52E- 


15 


6 


97E- 


11 


2 


82E- 


10 


3 


03E-11 


























5 


37E- 


12 




4p 


2 


25E- 


20 


3 


75E- 


16 


5 


73E- 


12 


2 


D3E- 


11 


2 


21E-12 


























5 


59E- 


13 




4d 


1 


15E- 


20 


1 


64E- 


16 


2 


40E- 


12 


8 


51E- 


12 


9 


35E-13 


































4f 


3 


41E- 


21 


7 


40E- 


17 


3 


93E- 


13 


7 


06E- 


13 


6 


59E-14 
































H 


-+Li+ 


3 


89E- 


12 


5 


37E- 


11 


7 


61E- 


08 


3 


47E- 


08 


3 


98E-09 


3 


30E- 


11 


7 


94E- 


12 





















Table 2. Predicted flux equivalent widths (in pm=10mA) for the 670.8 nm line in the Sun, HD140283 and HD84937 based on 
ID and 3D models, in both LTE and non-LTE. Results are shown for various model atoms as described in the text. Note that the 
3D results shown are based on only one snapshot and hence are not appropriate for a direct comparison with the ID cases. 

















ID 








3D 




Star 




Parameters 




log CLi 






^NLTE 


^NLTE 




^yNLTE 


^^NLTE 


^^NLTE 




Toff 


logs 


[Fe/H] 






noH 


wH 


wO . 5H 




noH 


wH 


wO . 5H 




[K] 








[pm] 


[pm] 


[pm] 


[pm] 


[pm] 


[pm] 


[pm] 


[pm] 


Sun 


5777 


4.44 


0.0 


1.1 


0.40 


0.34 


0.38 


0.37 


0.55 


0.37 


0.40 


0.39 


HD 140283 


5690 


3.87 


-2.5 


1.8 


2.40 


2.18 


2.66 


2.61 


3.84 


1.96 


2.35 


2.32 


HD84937 


6330 


4.04 


-2.25 


2.0 


1.31 


1.44 


1.57 


1.55 


1.79 


1.11 


1.26 


1.24 



cross-sections by 50% (wO . 5H). The remaining uncertainties 
in the charge transfer reactions are thus expected to have lim- 
ited impact on the Li I line formation. 

In terms of abundances, the inclusion of charge transfer 
reactions lowers the derived values by about 0.05-0. lOdex 
for all three stellar models, whether in ID or 3D. While 
the ID non-LTE abundance corrections (ID non-LTE - ID 
LTE) are small (<0.05dex) for the relevant stars (Table |2] 
see also Carlsson et al. .1994) . the 3D non-LTE effects (3D 
non-LTE - 3D LTE) are substantial for metal-poor stars 
(Asplund et al. '2003V However, since the differences between 
3D LTE and ID LTE are also significant and of opposite sign 



(Asplund et al.'1999), the net difference between 3D non-LTE 
and ID non-LTE is less dramatic. For the two metal-poor stars 
this difference amounts to about 0.05-0. lOdex lower abun- 
dance than in ID^. This has a direct bearing on the question 
of the primordial Li abundance. The observational evidence 
suggests a plateau at logeLi = 2.1 — 2.3 with ID non-LTE 

' Note that it is not appropriate to directly compare the equivalent 
widths in Table|2|to estimate the net abundance difference between ID 
and 3D as the 3D results are here only shown for one snapshot. The 
proper procedure is to add the 3D non-LTE corrections obtained here 
to the 3D LTE corrections published by Asplund et al. jl999t based 
on a long time-sequence with better numerical resolution. 
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-5 -3 -2 -10 1 2 

lOyioT 500nm 




-5 -4 -3 -2 - 1 2 

logiOT 500nrr, 



Fig. 1. Departure coefficients with standard optical depth for 
lower Li I levels in the solar ID model with log eu — 11 (up- 
per panel) and in the HD140283 ID model with log eu = 1.8 
(lower panel). Full lines are for the atomic model without H 
colUsions noH and the dashed lines are with H colUsions wH. 

abundance corrections taken into account (Ryan et al. 120001 
Bonifacio I2002t . with the abundance range partly reflecting 
different effective temperature scales for halo stars. By slightly 
decreasing these abundances, the inclusion of 3D models and 
charge transfer reactions slightly aggravates the discrepancy 
between the observed Li abundances and the predicted pri- 
mordial Li abundance from Big Bang nucleosynthesis. High- 
redshift deuterium measurements and analyses of the cosmic 
microwave background yield a relatively high baryon density 
of the Universe (e.g. Buries et al. l2001> . In particular, the re- 
cent WMAP analysis (Bennett et al. l2003> corresponds to a pri- 
mordial Li abundance of log eLi.p = 2.69 ±0.13 based on stan- 
dard Big Bang nucleosynthesis (Buries et al. .2001 ) . The small 
amount of scatter in observed halo Li abundances is very dif- 
ficult to reconcile with uniform stellar Li depletions of 0.4 dex 
or more. Likewise, it would require improbably high effective 
temperatures to resolve the disagreement. Clearly, halo stars 
and their spectrum formation are not as well understood as per- 
haps thought. 

4. Discussion 

We have found that while excitation and de-excitation by in- 
elastic H collisions are not important in thermalizing Li, the 
charge exchange reactions involving the Li(3s) state, are quite 



important. This mechanism is almost certainly important also 
in other alkalis (Na, K, Rb etc), where large cross-sections are 
also found (Dickinson et al. 1999; Janev & Radulovic lI978> . 
For other elements, further detailed calculations will be re- 
quired to determine whether inelastic H collisions, including 
excitation and charge exchange, are important. The main ob- 
stacle for such calculations at present is the quantum chem- 
ical potentials and couplings required for a large number of 
molecular states. Though systems of avoided ionic crossings 
exist in all metal hydride quasi-molecules, we note other mech- 
anisms may be important. As suggested by Lambert {1993}, 
the importance of details in the molecular dynamics, such as 
avoided ionic crossings, in this process makes a simple gen- 
eral recipe impossible, in contrast to collisional broadening 
and depolarization by H where such effects play a lesser role 
(Ans tee & 0'Mara .l995 : Barklem & 0'Mara l2Mn Derouich 
et al. l2003t . Ab initio potentials for alkaline-earth-hydrides are 
possible with current techniques (Chambaud & Lew Il989t . 
However, some cases of particular astrophysical interest such 
as O+H and Fe-nH present a considerable challenge. 
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